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Abstract 

New variables of separation for few integrable systems on the two-dimensional sphere with 
higher order integrals of motion are considered in detail. We explicitly describe canonical trans- 
formations of initial physical variables to the variables of separation and vice versa, calculate the 
corresponding quadratures and discuss some possible integrable deformations of initial systems . 



1 Introduction 

A fundamental requirement for new developments in mechanics is to unravel the geometry that un- 
derlies different dynamical systems, especially mechanical systems. There are several reasons why this 
geometrical understanding is fundamental. First, it is a key tool for reduction by symmetries and for 
the geometric characterization of the integrability and stability theories. Second, the effective use of 
numerical techniques is often based on the comprehension of the fundamental structures appearing in 
the dynamics of mechanical and control systems. In fact, geometric analysis of such systems reveals 
what they have in common and indicates the most suitable strategy to obtain and to analyze their 
solutions. 

Already in 19th century Euler and Lagrange established a mathematically satisfactory foundation 
of Newtonian mechanics. In |10j Jacobi united their ideas with the Hamilton optic theory and with 
the Abel geometric methods at a new Hamilton- Jacobi formalism. The Hamilton-Jacobi formalism 
was a crucial step towards Liouvilles classical definition of the notion of integrability [T5] based on the 
notion of first integrals of motion. 

The Liouville definition of integrable Hamiltonian systems naturally covered many classical ex- 
amples. Among them are the Kepler motion solved by Newton, harmonic oscillators solvable by 
trigonometric functions, the Euler and Lagrange spinning tops and the Jacobi example of geodesic 
motion on an ellipsoid solvable by elliptic functions [2J, the Neumann system on the sphere [16 and 
Kowalevski top |13j solved terms of hyperelliptic functions etc. In novel times much attention owing to 
the another discovery of the vast class of integrable soliton nonlinear partial differential equation, that 
admits this type of integrability when dynamics is restricted to finite dimensional Liouville tori and 
the system appeared to be completely integrable in the Liouville- Arnold sense. They all are more or 
less connected with the hyperelliptic curves and with the hyperelliptic functions [TJ [TUJ [T7] . Below we 
show that foregoing development of the theory detected a number of cases when associated algebraic 
curve is non hyperelliptic and and its genus exceeds the number degrees of freedom [24l [26j [27] . 

Bi-Hamiltonian structures can be seen as a dual formulation of integrability and separability, in 
the sense that they substitute a hierarchy of compatible Poisson structures to the hierarchy of functions 
in involution, which may be treated either as integrals of motion or as variables of separation for some 
dynamical system [2"0] . The Eisenhart-Benenti theory was embedded into the bi-Hamiltonian set-up 
using the lifting of the conformal Killing tensor that lies at the heart of Benenti's construction, which 
may be realized as a computer algorithm [9]. The concept of natural Poisson bivectors allows us to 
generalize this construction and to study systems with quadratic and higher order integrals of motion 
in framework of a single theory [20, 23, 24 . 
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The aim of this note is to discuss separation of variables for integrable natural systems on the 
two-dimensional unit sphere § 2 from [31] [33J 1231 (Ml HI] • In the above mentioned previous papers we 
focused our attention on the bi-Hamiltonian calculations of the variables of separation starting from the 
given integrals of motion. This note is devoted to construction of the initial physical variables in terms 
of variables of separation, to calculation of the corresponding quadratures and to discussion of the 
possible integrable "gyroscopic" deformations of these systems associated with genus three algebraic 
curves. 

In order to describe integrable systems on the sphere we will use the angular momentum vector 
J = (Ji, J2, J3) and the Poisson vector x — (x%, x%, X3) in a moving frame of coordinates attached to 
the principal axes of inertia. The Poisson brackets between these variables 

{ Ji ) Jj 1 — &ijk Jk j { Ji ) •Ej } — ^ijk-^k ) {^i i } — j (I'l) 

may be associated to the Lie- Poisson brackets on the algebra e*(3). Using the Hamilton function H 
and the Lie-Poisson bracket {.,.} (|1.1[) on the Euclidean algebra e*(3) the customary Euler-Poisson 
equations may be rewritten in the Hamiltonian form 

ji = {Ji,H}, Xi = {x h H}. (1.2) 

Remind, that the Lie-Poisson dynamics on e*(3) can be interpreted as resulting from reduction by the 
symmetry Euclidean group E(3) of the full dynamics on the twelve-dimensional phase space T*E(2>) 
[3J. There are two Casimir elements 

3 3 

d = |af = 5> 2 fe , C 2 Hx,J) = 5> fe J fe , (1-3) 

fc=l k=l 

where (., .) means inner product. Using canonical transformations x — > ax we will always put C\ = 1 
without loss of generality. 

If the square integral of motion C2 = (x, J) is equal to zero, rigid body dynamics may be restricted 
on the unit sphere S 2 and we can use standard spherical coordinate system on its cotangent bundle 
T*S 2 



xi = sin sins' , X2 — costpsmt', X3 = cos( 

sin <j> cos 9 cos cf> cos 8 

Jx = ^-p. — p<t>- cos (pp e , J 2 = r-7 — p0+sm0p0, J3 = -p<f> 

smO smO 



(1.4) 



We use these variables in order to determine canonical variables of separation on T*S 2 . 

As usual all the results are presented up to the linear canonical transformations, which consist of 
rotations 

x^>aUx, J —> U J , (1.5) 
where a is an arbitrary parameter and U is an orthogonal constant matrix, and shifts 

x — y x , J— >J + Sx, (1-6) 

where S is an arbitrary 3x3 skew-symmetric constant matrix. 

Of course, any canonical transformation of the spherical variables (ll.4[) yields automorphism of 
e*(3) too. For instance, trivial canonical transformation 

P6^Pe + f(0) (1.7) 
gives rise to "generalized" shift depending on arbitrary function /(X3): 

Jt-ji-^p^, w^J^SsL, (i.8) 

\JX{ + X2 \Jx{ + X2 

This and more complicated canonical transformations of e* (3) are discussed in [2j [11] . 
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2 Kowalevski top and Chaplygin system 

Following to [HI [HI [53], we determine canonical coordinates qi j2 on T*§ 2 as roots of the following 
polynomial 

pi sin 2 9 + p 2 , cos 2 9 „ 
B(A) = (A- gl )(A-g 2 ) = A 2 - ^ , ^ A-« 2 -fr 2 

sin 6 1 cos^ 9 

(2.1) 

(a cos ad> - b sin a</>) (p 2 , sin 2 + p 2 cos 2 6) 2 sin 6»(a sin ad> + b cos a^)p^ Pe 
sin™ # cos 2 sin" 9 cos 2 # 

Then we can introduce auxiliary polynomial 

sin 9(a cos ad) — b sin ad) 



A(X) = ^A + ^— 21 P0 

a cos £/ a 



such that 



{B(\),A0*)} = j^— ((/i 2 -a 2 - b 2 )B{\) - (A 2 - a 2 - b 2 )B(^)) , {A(A), A( M )} = . 



It entails that 

Pj = ~ u 2 A ( X = ^' J' = 1 ' 2 ' ( 2 - 2 ) 

are canonically conjugated momenta on T*S 2 with the standard Poisson brackets 

{Qi>Pj} = hj , {iuQa} = {puPi} = o. 

Below we prove that at a = 1, 2 this variables are variables of separation for the Kowalevski top and 
Chaplygin system, respectively. 

2.1 Kowalevski top 

Let us consider Kowalevski top defined by the following integrals of motion 

Hi = J 2 + J 2 + 2 J 2 + 2bxi (2.3) 



H 2 = (J 2 + J 2 ) 2 - 4b(xi{J 2 - J 2 ) + 2x 2 JiJ 2 ) - 4b 2 x 
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In original Kowalevski work [T3] the first step in the separation of variables method is the complexifi- 
cation: she introduces 

zi = Ji + iJ 2 , z 2 = Ji — iJ 2 
as independent complex variables. Next she makes her famous change of variables 



_ fl(zi,Z 2 ) ± V /fl(zi,Zl)fi(z2,Z 2 ) 

51,2 " 2( Zl -z 2 ) 2 
The fourth degree polynomials R(zi, z^) may be found in [13l[2]. It brings the system (| 1 . 2|) to the form 

(-l) fe (si - s 2 )s k = VP(st) , fc = l,2, 

where 



P(s) =A((s-Hi) 



i2 



H 2 + 4b 2 Ci 



s[{s -Hi) 2 + b 2 Ci- H2 + * a Cl )+b 2 C 2 



(2.4) 
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Consequently, initial equations of motion can be written as hyperelliptic quadratures 



glfl , S 2 S 2 

where we can substitute the conjugated momenta p Sk instead of y/P(sk) in order to get standard 
Abel-Jacobi form. So, the problem can be integrated in term of genus two hyperelliptic functions of 
time. Finally, we have to substitute these functions of time Sfc(t) and p 8k {t) into the initial variables 
x, J, the corresponding expressions may be found in [131 112] . 

Discussion of the another variables of separation for some particular subcases in the Kowalevski 
dynamic may be found in [2 . As usual different variables of separation are related with the distinct 
integrable deformations of the initial integrals of motion. 



2.1.1 New real variables of separation at C% = 

According to [23], at a = 1 coordinates q± i2 (|2.ip are variables of separation associated with the 
Hamilton function 

H = Jl + Jf + 2J| + +2ax 2 + 2bxi , 

which may be reduced to the initial Hamiltonian Hi using rotations (jl.5|) around the third axis 
so we can put a = in (|2.1[) without loss of generality. 

Coordinates qx,i (EHJ at a = 1 and a = are defined by 

12 T 2 I j2 , j2\ b(xi(J 2 -J 2 2 )+2x 2 J 1 J 2 ) 

B(X) = (A — 9l )(A - ga) = A 2 + + X ^ 2 Jl + - h) -A J 

The conjugated momenta are equal to 



A(A = q%) , . x x J 2 -x 2 J 1 + x\ J 2 

P k = 2 — rr~ ' A ( A > = A H • 

q% - o £3 x 3 

This variables are differed on a constant terms ±6 from variables introduced in [22 . Inverse 
transformation reads as 

„2^2 , IU 2 „2^„2^ 2( b 2 - ql)(b 2 - q 2 ) 

Km - 12) 2 



— 91^2 /,, 2 2\ 2 1 /.2 2\ 2 . 

Xl = bJqT^^V -1l)Pl + ( b -1 2 )P 2 ) — ~ 1>U>2 



x 2 



x 2 [ b (pi -P2 2 - (pi9i -p 2 g 2 

%l - 92j 2 V 

r (^_ g 2 )2p 4 + (6 2_ g 2 )2p 4 ^ 2(6 2 -g 2 )(6 2 -g 2 )p"frl 

V (91 - ?2) 2 (<?i - ?2) 2 



(2.5) 



Ji 



Jo = 



(6 2 - q\)p\ - {b 2 - qfjpl b 
92 Q 2 ~ q 2 )pi ~ qi{b 2 - q 2 )p 2 X3 

(6 2 _ ^2 _ (6 2 _ g 2 )p 2 ft . 



J 3 = -\{ql-V){b 2 -ql) 



qi - 92 

Coordinates of separation take values only in the following intervals 

qi > b > q 2 , 
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similar to the standard elliptic coordinates on the sphere [2]. 
In this variables integrals of motion H\ 2 (|2.3p look like 



Hi 



H 2 = 



(b 2 -q 2 ) 2 pi-(b 2 -q 2 ) 2 pj-(q 2 -q 2 ) 
(b 2 - q \) V \ - (b 2 - q\)v\ 

({b 2 - ql)vl - (b 2 - q\)p\ + 9i + 92) ((& 2 - q\)v\ ~ (b 2 ~ q\)vl - qi - 92) 
(b 2 - q\)v\ - (b 2 - qWl 

x ((6 2 - q\)p\ - (b 2 - ql)p% + 9l - q 2 ) ({b 2 - q\) V \ - (b 2 - q 2 2 )p 2 2 - q x + q 2 ) . 

It is easy to see that integrals of motion and variables of separation are related via the following 
separated relations 



<I> 



(2{q 2 - b 2 )p 2 +H X + v/fla) (2(g 2 - b 2 )p 2 + H X - y/H^j - M 2 = 0, (2.6) 



at q = qi, 2 and p = pi,2- Equation $(q,p) = defines genus three hyperelliptic curve with the following 
base of the holomorphic differentials 

n x = * ri a qdq 



do 



p(b 2 - q 2 ) (Hi -2(b 2 - q 2 )p 2 ^ ' p(b 2 - q 2 ) [h x - 2(b 2 - q 2 )p 2 

Pdq 

ffi - 2(6 2 - g 2 )p 2 ' ' 



In fact equation (|2.6[) is invariant with respect to involution (q,p) — > (—q,p). Factorization with respect 
to this involution give rise to elliptic curve. 

In variables of separation equations of motion (|1.2[) have the following form 



qi qi 



Pl (b 2 - q 2 )^, 2(b 2 q 2 )p\) p 2 {b 2 - <&)(Ht 2{b 2 q 2 )p 2 



= 0. 



9i , 92 _ 2 _ 



H 1 -2(b 2 -q 2 )p 2 H 1 -2(b 2 -q 2 )p 2 

The above quadratures in the integral form 

gi rii rqi rqi 

Qi = fa , / fl 3 + n 3 = -2t + f3 2 , (2.7) 
go Jqo J qa J qo 

represent the Abel-Jacobi map associated to the genus three hyperelliptic curve defined by &(q,p) = 0. 
In particular it means that instead of p in (|2.7p we have to substitute function on q obtained from 
the separated relation (|2.6[) . 

In order to give explicit theta-functions solution one can apply some remarkable relations between 
roots of certain functions on symmetric products of such curves and quotients of theta-functions with 
half-integer characteristics, which are historically referred to as root function and are generalized so- 
called Wurzelfunktionen that were used by Jacobi for the case of ordinary hyperelliptic Jacobians 
P2 [6] . For the case of odd order hyperelliptic curves such functions were obtained by Weierstrass [28] . 
Inverting the map (|2.7[) and substituting symmetric functions of q\,q 2 ,p\,p 2 into (|2.5[) . one finally 
finds x, J as functions of time. 

2.1.2 Deformations of the Kowalewski top 

According to [2H [H] , using separated relations 



$1 = [2{q 2 -a 2 )p 2 + H 1 + yH 2 )[2(q 2 -a z )p 2 + H 1 -yH 2 ) - 4cu z + Adu + e(q 2 - a 2 )p = , (2 
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one gets Hamilton function of the generalized Kowalevski top 

- C -^) (J! + 4) + 2JI + 2t Xl + - 7 =1= + f - gl / a) 

Second integral of motion is equal to 



fr (ttj+C-l) 2 2 2 2 / 4(X! - H J| + 2X 2 J X J 2 )(x 2 + C - 1) ^^1+^1 

#2 = -4 (J1+J2) - I 



4(^ + c- l)a 2 



a?3 x 3 

' (2Q^ + J g)(c-£^--£|) (x 2 Jx -x 1 J 2 )e \ rf 2 

v/^f+ x 2 a: 2 + 2(x 2 +x 2 )x 3 J + x\ + x\ 



( Ji + Jj)(x2Ji - xi J 2 )(c - a; 2 - ^)e (x 2 Ji - XiJ 2 ) 2 e 2 

2 v ^f+^|x3 3 16(x 2 + x 2 )x 2 1 ' ' 

According to [TTJ HI] canonical transformation (|1.8[) reduces Hamilton function (12.9[) to the natural 
form 

M 1 " ^ ) (J ' + JJ) + 2J| + 2 "' + vsfcf - figFR ' (2 - u) 

at 

/(X3) = 8(, 2 + c-l)- 

At c = 1 this system coincides with one of the deformations discussed in . Below we will show only 
the final form (|2.11[) of the deformed Hamiltonians and will omit the intermediate form (|2.9p for the 
brevity. 

It is easy to calculate the corresponding equations of motion 

eh h __ 



(b 2 - q 2 )[8Hx P x + e - 16p?(fe2 _ 5 2)j p2(6 2 _ q j)(8H lP2 + e - lGp 3 ^ 2 _ ^ 

9i ^ 92 = _1 

SHipi + e - 16p?(6 2 - gf ) 8H lP2 + e - l6 P 3 2 (b 2 - q 2 ) 4 

and prove that the Abel-Jacobi map on genus three hyperelliptic curve has the same form (|2.7[) 
qi rqi rqi rqi 

sii+ tlx = pi, / n 3 + / n 3 = -2t + (3 2 , 



iq Jq Jqo J qo 

where P have to be solution of the separated relation (I2.8P and 



{b 2 - q 2 )(8H lP + e - l6 P 3 (b 2 - q 2 )) 8H lP + e- l6 P 3 (b 2 - q 2 ) 

2.2 Chaplygin system 

Let us consider Chaplygin system defined by the following Hamilton function 

Rx = J? + Jf + 2J| - 2a{x\ - x\) - 2bxxx 2 - . (2.12) 

^3 

At c = this system and the corresponding variables of separation have been investigated by Chaplygin 
[JJ . Singular term has been added by Goryachev in [8] . 

Using rotations (|1.5p around the third axis [TT] we can put 6 = without loss of generality. In 
this case second integral of motion is equal to 



c 



2 



H 2 = Jf + Ji- — ) - 4ax 2 (J 2 - J 2 ) + Aa 2 x 4 3 . 
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According to [2lJ 121], coordinates gi,2 (|2.ip are variables of separation for this integrable system at 
a — 2 and b — 0. In this case gi l2 are roots of the following polynomial (|2.ip 

t2 , 7 2 9o 7 2 or T 2 4- 7 2 1 

x 3 x 3 ^ 

whereas momenta px 2 are values of the other auxiliary polynomial 

_ x 2 Ji - xiJ 2 . _ C1X1X2J3 _ a(x 2 - ggXxgJi " J 2 ) 
2a; 3 x\+xl 2{x\ + x%)x 3 

at A = (71.2 (|2.2j) . Inverse transformation reads as 



Xi 



J'2 



v /2(g 1 - a)(q - g 2 )(pi(gi + a) ~P2(g2 + a)) 
\/a(oi - 92) 



y/2(q 1 + a)(g 2 + a)(pi(gi - a) -^2(92 - «)) 
Va(qi - 02) 



„ _ t - 4 M --'»"' -M-ra (2 . 13) 

V 01 - ?2 



(a + gi)(a + g 2 ) , {qi~a){a-q 2 ) 
Jl = ^ 2^ X3 ' J2 = -V 2a" " 3 



J 3 = -2 J (g 2 - a 2 ) (a 2 - g|) — — — 
v gi - g 2 

As usual coordinates of separation take values only in the following intervals 

gi > a > g 2 • 

This variables gi,2 are related with variables of separation from |21] by the rule qk — > qk + a,. 
In variables of separation integrals of motion read as 

H l = 4 (« 2 - 1l)Pl + 4 (« 2 - iDvl + H + 12 - -T7— 2 20 yA 2T^~ ' 

4(a 2 - q{)p{ - 4(a 2 - q 2 )p 2 + qi - q 2 



(4(a 2 - ql)p\ - 4(a 2 - q\)p\ + q x - g 2 ) + 

[A{a 2 - q 2 )p 2 - 4(a 2 - gf> 2 + q x - g 2 J 
- 2c(gi + g 2 ) . 

It is easy to see that integrals of motion and variables of separation are related via the following 
separated relations 

,2\„2 



$ = 8(q 2 - a 2 )p 2 -2q + H 1 - ^ H 2 8(g 2 ~ a 2 )p 2 -2q + H 1 + ^ H 2 - 4cg = 0, (2.14) 



at g = gi,2 and p — p\. 2 . Equation $(g,p) = defines genus two algebraic curve with the following 
holomorphic differentials 

da U{a 2 -q 2 )p 2 + q\dq 

fli= 7 V' n2 = —/ " V ( 2 - 15 ) 

p(a 2 - q 2 ) (Hi - 8(a 2 - g 2 )p 2 - 2gJ p{a 2 - q 2 ) ( H x - 8(a 2 - g 2 )p 2 - 2gJ 
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The corresponding quadratures look like 

9i <72 



Pl (a 2 - qf^Hx - 8(a 2 - q\)p\ - 2q x ) p 2 (a 2 - q*)^ - 8(a 2 - q 2 )p 2 - 2q 2 
(4(a 2 - q\) V \ + qi) qi (4(a 2 - q 2 )p 2 + g 2 ) <?2 



Pl (a 2 - q 2 )^ - 8(a 2 - q 2 )p 2 - 2 Ql ) p 2 (a 2 - q^Hx - 8(a 2 - q 2 )p 2 - 2q, 
The Abel-Jacobi map on genus two hyperelliptic curve has the standard form 

fQ2 rqi rq2 

fti + / Sli = Pi, / ^2 + / n 2 =8t + (3 2 , 

q Jq Jqo J q 

where p into fii 2 have to be solution of the separated relation (12 . 14[> . 
2.2.1 Deformations of the Chaplygin system 

According to [2lJ[24], if we substitute this variables of separation onto the following separated relations 



$1 = (d,(q 2 -a 2 )p 2 -2dq + H 1 -\l H 2 ^(8(q 2 -a 2 )p 2 -2dq + H 1 + y H 2 ^-4cq + e{q 2 -a 2 )p = Q, (2.16) 
one gets the Hamilton function of the generalized Chaplygin system 

ft - ft - ^) + ■*> + - MA - 4) - »«, - j-^ + ■ (2 ' 17) 



As for the Kowalevski top, using canonical transformation (|1.8I) at 

/(*a) 



16(d- l + x 2 ) 2 

we can reduce Hamilton function (|2.17l) to the natural Hamiltonian 

Hi = ( 1 - (J? + Jf ) + 2 J 3 2 - 2a(a; 2 - xl) - 2bx lX2 ' ' ' ' 



. 2 ; VJlT,J2 ^^ 3_ ^ 1_ - i2 ^^ i - i2_ rf-l + a ; 2T 256(d-l + a ; 2 )3- 

At d = 1 additional term is equal to e(a;^" 4 — x^ 6 ) and this system coincides with one of the deformations 
considered in [21?] . 

In this case we have genus three hyperelliptic curve with holomorphic differentials 



(a 2 - 9 2 ) (e + 32p(#i - 8(a 2 - <? 2 )p 2 - 2dq)} 

qdq 

(a 2 - q 2 ) (e + 32p(7Ji - 8(a 2 - q 2 )p 2 - 2dq) 

p 2 dq 



e + 32p(H 1 - 8(a 2 - q 2 )p 2 - 2dq) 
and the corresponding quadratures involve all this differentials 

rll r<l2 rqi rqi j- 

/ n x + / rti = pi, / {4n 2 + dn 3 )+ (4f> 2 + df> 3 ) = -- + 

•/go •'90 •'go •'go 

in contrast with other integrable systems on genus three algebraic curves considered in this note. 



3 Integrable systems associated with trigonal curves 

According to [251 [Ml HZ], we introduce other coordinates qi j2 on T 8 § 2 defined as roots of the following 
polynomial 

B(X) = (X-qi)(X-q 2 ) = A 2 - i^FA + A , i = V=I, (3.1) 



with coefficients 

F = (g{6) Pe - ih(6)p^j * , A = a exp U<f> - J d6) , (3.2) 
depending on arbitrary functions g(9) and h(6). As usual conjugated momenta p\^ 2 are equal to 

p k =A(X = q k ), A (X) = ij^- ) - i f. (3.3) 
It is easy to prove, that these polynomials satisfy to the following relations 

{B(X), A(m)} = ~~~r — ) , {A(X), A(m)} = {5(A), = , (3.4) 



/i — A \ A /i 
which give rise to canonical Poisson brackets 

Ui,Pj} = kj , {qi;Q2} = {Pl,P2} = 0. 

Substituting variables 

x = aq k 1 , z = a p kl k = l,2, a, a GlR, (3.5) 

into the generic equation of the (3,4) algebraic curve 

$(z, x) = z 3 + {aix + a 2 )z 2 + {H x x 2 + b\x + b 2 )z + x* + H 2 x 3 + c x x 2 + c 2 x + c 3 = , (3.6) 

and solving the resulting equations with respect to Hi 2 , one gets the following Hamilton function 

/c 2 + ia 6iw 2 - ao a i w 2 . , 2a aiw 2 - i&i A . gw 2 {c 2 + ia Q b 1 w 2 - a^aiiwf) 

ti\= 1 + V + a H ip0 , 

\ aoaw 2 aw 2 J aao 

(3.7) 

where geodesic Hamiltonian T and potential V are equal to 

a o(^ 2w 2 — 3/w;2 + 3) ia a 2 (hw 2 — l) 2 b 2 h(hw 2 — 1) ic 3 /i 2 \ 2 
a 2 a 2 w 2 a 2 W2 a^a 2 w 2 J ^ 

i-9 Z' fo_2...3 o:_ _ ...2 01. 1 2ic 3 



/ 2 3 2 ^l c 3 \ 2 2 \ 

\2a w 2 — 2iaoa 2 w 2 — 2b 2 w2 H J /i — 3a u; 2 + 2ia a 2 W2 + b 2 J p^Pe 



a 2 w 2 \\ - ~ a 

g 2 (a b 2 w 2 + iala 2 wl - a^w 3 . - ic 3 ) 2 

+ 9 

a z aow 2 

ia 2 e~'^ ( (ao&2W 2 + iaga 2 W2 — a^w 3 , — ic^awie 1 ^ | ici 



Here 



aaowiu'2 aoa 2 w 2 qqw 2 



wi = exp — / «o , w 2 — ' 



Second integral of motion H 2 is a cubic polynomial in momenta p^ and pg. 
The resulting Hamiltonian H\ (|3. 7|) has the natural form, if and only if 

2a a 1 w 2 — i&i = , c 2 + ia 6iw 2 — apa!?^ = . 



So, because W2 ^ 0, we have to put 

a± = h = c 2 = . 

If we want to obtain diagonal metric, then we have to solve integral equation 
2h(aQW% — \a^a,2Vj\ — 00^2^2 + IC3) — Sa^w^ + 2iaga2W2 + dob- 



0. 



(3.8) 



with respect to functions h(9),W2(0) and parameters 00,(22,02,03. If we want to get real potential 

V = /i(0)cos($ + / 2 (0) 
in p.7p . we have to add one more equation to (|3.8|) 

ia 2 (ag?i'2 — ia^a^u^ — ao02U>2 + ic3)u> 2 + a 4 = (3-9) 



depending in addition on function w\ and parameters a (|3 . 5|) and a (I3.2[) . 

Some particular solutions of these equation have been studied in [THl [23 121] including integrable 
systems due to Goryachev, Chaplygin, Dullin, Matveev etc. For all these systems, we collect a and 
the zero-valued coefficients in (|3.8[) in the following table 



Goryachev-Chaplygin top 


ao 


= 2ia 


62 


= c 3 


= 


Goryachev system 


a 


= 2ia/3 


02 


= °2 


= 


Case 3 from [TH] 


a Q 


= ia/3 


a 2 


= b 2 


= 


Dullin-Matveev system 


a Q 


= ia 


C3 


= 




Case 5 from [19] 


a 


= ia/2 


a2 


= c 3 


= 



Integrable systems with the same coefficients in the separated relations (|3.6[) and with different ao and 
a' Q in (13.5[) are related by non-canonical transformation of the momenta 



z = a Q p k -> z = a p k , 



k= 1,2. 



(3.10) 



3.1 Goryachev-Chaplygin top 

Let us consider Gory chev- Chaplygin top with the following integrals of motion 



Hi = J'l + Jf + 4J3 + azi 



#2 = 2 J 3 ( Jl + J| 



0x3 Ji 



(3.H) 



In this case variables of separation (I3.1I3.3|) are determined by 

Ji + i J 2 



qi + <?2 



2J3 

x\ x 3 (xi + IX2) 



HQ2 



2x\{xi + ix 2 ) 



Pi, 2 



91,2 
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They are related with initial variables by the rule 



ia(q 1 -q 2 ) / 2i(piqi - q 2 p 2 ) 

Xi+lX-2 = — -. r, X 3 - 



Jl + 1J2 



X\ - IX 2 



49i9i(pi9i - 92P2) ' V <h - 92 

a(<7iPi - t mq2(pi -pi) 



2qiq 2 {Piqi - 92P2J1 



9l - 92 

49192 



«(9i - 92) 
a(9i - Q2) 2 \l- 



> ((i " 2pi)ofpi + (4pip 2 - ipi - 1^2)9192 + (i - 2p 2 )q\p 2 ^ 



T -j ^9192 / 2 

Jl-lJa = , =7 F= I glPl ~ g2P2j(l- 2^1)^X9! + 

/ 2i(pigi - g 2 P2j v 



9i - 92 

+ (9iPi + 92P2)(i - 2p 2 )p 2 ql 
Separated relation is given by equation with real coefficients 

4 

Equations of motion in variables of separation look like 

9i 92 



A*) = (p 2 -b)q 2 + (p 3 - H 1 fi + H 2 )q+ i- =0, 9 = 91,2, /i = 2i 91,2^1,2 ■ (3.12) 



= 0. 



2i. 



9i(3/z 2 + - Fi) 92(3^i + 2g 2 M2 - #1) 

Mi9i ! ^292 

91 (3/if + 291/ii - ffx) g 2 (3ju| + 292M2 - #1) 

By making the birational change 

q = h » = l (3 - 13) 

the curve (|3.12[) can be transformed to the canonical trigonal form (|3.6p at 

ai = 61 = b 2 = c 2 = c 3 = , 
whereas other parameters are functions on a, b. 

3.1.1 Deformation of the Goryachev-Chaplygin top 

Substituting q = qi 2 and /1 = 2igi i 2Pi 1 2 into the non-hyperclliptic algebraic curve of genus three 
defined by the following equation 

a 2 

(i) = eg 3 + (/i 2 + d[i-b)q 2 + (p 3 + ep 2 - H lf i + H 2 )q + — = 0, (3.14) 

and solving a pair of the resulting equations with respect to Hi 2 one gets deformation of the initial 
Hamilton function 

~ t2t2.t2 b ( c—d + e c 2c \ T (c-dxl+exi) 2 



"3 

using the generalized shift f| 1 .8[) at 



/ = - 



-r>2 I t»2 ™ 2 ryi4 / /| ™U / ,--vi 2 I ,71 2 \ 

•j. ^ 1^ iV 2 3 3 / 3 V 1 2/ 



i(ex| — dx§ + c) 



/yj 2 /y* 3 

0,3 J, 3 
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In this case quadratures arc defined by the following differential equations 

2 

— ; — 5 — 7T— = , fJ.k=2i q k Pk 

k= i qk{i^l + %Pkqk + 2e/i fe + dq k - Hi) 

V ^ ^- = 2i. 

k=1 Qk(3lJ>l + Z^kqk + 2efik + dq k - Hi) 

If c = and d = e, one gets the usual Goryachev-Chaplygin gyrostat with the Hamiltonian 

Hx = J\ + Jf + 4 J 3 2 - e J 3 + axi + -g- ■ 
In this case equation (|3.14p defines genus two hyperelliptic curve instead of trigonal one. 
3.2 Goryachev system 

Let us consider Gorychev system on the sphere defined by the following integrals of motion 

2 2^2 ax i b 
Hi = J 1 + J 2 + — J3 1 



Q & ' 2/3 2/3 ' 

V.J T.J 



tt 2 J 3 / 2 2 8 2 6 \ a(3x 3 Ji - 2a;i J 3 ) 

= "~3"V 5 5F; — MP — ' 

The corresponding variables of separation q\.2 and pi : 2 (|3.1|3.3|) are obtained from 

4/3 T -i T T n 1/3 Q . 2/3 . 7 

qi+q2 = y L 2+ \ 2 ~ ' <7i<72 = ^7 ~- — V' Pl,2 = + 

1 — 2:3 1 - x 3 2(xi+ia;2j 2 g li2 

Inverse transformation looks like 

a 7 i<7i<72(pi -Vi) 

xi + ix 2 = , J3 = 



2<?i<3 , 2 ' qi - 32 

gig 2 (l-a;§) . + 92) 

Xl — 1^2 = 2 , Ji+lJ 2 = jV, 

a 2qiq 2 x 3 ' 

T -j ^hl<ll{Pi ~ P2) . 2qiq 2 {l - xl)(qi + q 2 ) 

Jl — 1J 2 — z r 131 775 

a (9i - 92) ax 3 f 
2i(pxqi -p 2 <72) X 

X3 ~ 



where 

/ 2i 

.Tt = I 

3(5i - q2 

Separated relation is given by equation with real coefficients 

a 2 2i 
®{q, M) = 9 4 - fo 9 2 + (M 3 - #1^ + H 2)q + — = , at g = g 1)2 , ,u = — qi, 2 Pi,2 

In this case quadratures read as 

rqi dq f g2 dq 



„„ <z(3 M 2 -#i) J qo q(3p 2 -Hi) 

pdq f q ' 2 pdq 2i 



2 ■ 



g(3 M 2 -ffi) 7 9Q q(3p 2 -Hi) 3 
As usual, here /i is a function on g obtained from the separated relation (|3.17l) . 



12 



3.2.1 Deformation of the Goryachev system 

Using trigonal curve of genus three defined by the following equation 

,2 



$1 = q 4 + cq 3 - bq 2 + (p 3 + dfl 2 - H x p + H 2 )q + = , (3.19) 



instead of (|3.17p one gets deformation of the initial Hamilton function p. 151) 

d | d + cx 2 / 3 \ j ^ (c + dx^ 3 ) 2 

xf+xl J _|_ ^2^2 



after the generalized shift (|1.8[) at 

i(c + dx^ 3 ) 



/ = - 



2 I/ 3 



3^3 



The corresponding equations of motion look like 



9i 92 „ 2i 



gi (3/x? + 2d Ml - ff- ) q 2 (3^ 2 + 2dfi 2 -H 1 ) 3 

/fi9i ^ Mi 92 _ 2i 

gi(3/xf + 2d/ii - Hi) q 2 (3f4 + 2d/i 2 - Hi) 3 ' 

Birational transformation (|3.13[) maps the curve (|3.19p to the canonical trigonal form (|3.6p at a 2 
b 1 = b 2 = 0. 

3.3 Case 3 from [19] 

Let us consider one more integrable system from [H]] defined by the following integrals of motion 
Hi = jf + jf+fl + g^Ei^J? ' ^ . b 



12 2(x 3 + l)J 3 (.T3 + I) 5 / 6 (xa + l) 1 ^' 

(3.21) 



H, = ±- J| - I JsH! - a(x 3 + 1)V8 J, + " ! /; 



27 3 3 "- 3 " 1 ' ~' ~ 1 ' 2(^3 + l) 5 /6 

The corresponding variables of separation 512 and pi t2 (|3 . 1 13.3[) are obtained from 

(l + a: 3 ) 2/3 J 3 i(x 2 Ji - xi J 2 ) a^/l + x 3 1/3 iJ 3 

9i+92 = 7T, 7\ 77- 7773-7 77. 9i92 = -^ —. — r, Pi, 2 = 31(1 + 13 )' + 

2(^3-1) (1 + xzyi A {xj, - 1) 2(asi+ia;2) 91,2 

Inverse transformation looks like 



ay/1 + x 3 19192 (Pi -P2) 

xi + 1x2 = — , J 3 - 



29i92 91 - 92 

29192(2:3-1) t 1*7 ia(pi-p 2 ) a(gi+ 92 ) 

zi - 1x2 = 7= , Ji + 1J2 = 



ay/TTx^ ' 4(gi - g 2 )\/lT^3~ 2< Zl g 2 (l + X3) 1/6 ' 

j • t _ 19191(3x3 + l)(pi -p 2 ) 2gig 2 (gi + 92) (x 3 - 1) 99 ^ 

1 2 " aCft-ftJVrr^ a(l+x 3 )V6 ' 

where 

_ i(pi9i ~P2g2) 3 _ 
X3 " 27(g 1 -< ?2 ) 3 
Separated relations are defined by equation with the real coefficients 

2 

<5>{q,lJL) = 2 q A -bq 2 + {^q-H 1 n + H 2 )q+ ±- = , q = g 1>2) fi= 1 ^^ ( 3 . 23 ) 
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The corresponding quadratures are given by 



= ft 



qi^-Hx) J qo q$f-Hi) 

(3.24) 

m , f q ' 2 m i . , a 



<lu 



q(3^-H 1 ) J qo q{ZlL 2 -Hi) 3 



As for the Goryachev system, birational change (I3.13|) transforms the equation (I3.23[) to the canonical 
trigonal form (13.6[) at 

ai = «2 = bi = b 2 = c 2 = , 

It allows us to prove that integrals of motion for this system (13.211) are related with integrals of motion 
()3.15|) for the Goryachev system by the non-canonical transformation (13.10P . 

It may seem that quadratures (|3~T8]) and (f3T24j) are trivially related by change of time 

but we have to keep firmly in mind that /i in (|3.18p is a function on q obtained from (|3.17p . whereas 
fj, in (|3.24p is another function on q obtained from (|3 .23[) . 

3.3.1 Deformation of the system (|3.21|) 

Similar to the Goryachev system, we can add two term to the initial trigonal curve of genus three 
(l3~23l) 

a 2 

$1 = 2q 4 + cq 3 - bq 2 + + d^ - H lf i + H 2 )q + — = . (3.25) 
Deformation of the initial Hamilton function (I3.2ip looks like 

Si = Hi (*__*__ cd W 3 ) J3 + (<vrr^+c(i + x 3 )-v«) a (3 _ 26) 



6 x 3 -l 2(a 3 -l) / 4(1 -a? 3 ) 

after canonical transformation (jl.8l) at 

i(d(l+x 3 ) + c(l + a: 3 ) 1 /3) 



/ = - 



In this case equations of motion are equal to 



2 V /I^f 



, = 7; QkPk 



qi (3^? + 2dfn - Hi ) q 2 (3yu| + 2dfi 2 - fli) 3 

Migi ! /i292 _ }_ 

qi(3f4 + 2dfi 1 -H l ) q 2 (3^ + 2dfi 2 - H 1 ) 3 ' 

3.4 Dullin-Matveev system 

Let us consider the Dullin-Matveev system [5] defined by the following integrals of motion 

X3 xi — \x\ 2 \ 2 ax\ 



(3.27) 



X 3 + C 4(X 3 + C) 2 J {x^ + c) 1 / 2 



(3.28) 



2 , 7 2 J3 , (4x 2 + 6x 3 c + c 2 + \x\ 2 )J 2 



H 2 = -{Jt + J: 
+ a V x 3 + c Ji 



2 4 4(a; 3 + c) 2 a; 3 + c 

axi J3 
2V^3 + c 
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According to [24] variables of separation 51,2 and pi )2 are defined by (|3.f I3.3[) 

J 3 Ji +ij 2 a iJ 3 

91+92 = -^- r -: — , 9i92 = ^7 -; — w , Pi, 2 = i(c + X3) H 

or by inverse transformation 

a i(pi9i - 92^2) t igi^Cpi -P2) 

Xi + 1X 2 = ; = . , X 3 = C , J 3 = , 

/ i[piqi -piqi) 9i - 92 91-92 

2W gig 2 

V 9i-92 



2\/ 9i92 



X-L-1X2 = -p- ( (c+ 1 + ipi)qi - (c+ f + ip 2 )92 

a(9i - 92)^ 



((c - 1 + \pi)qi - (c - 1 + ip 2 )g 2 J 



T , . T a (9i (29i + ggjgl - 92(2g 2 + qi)p 2 ) ,„ „„, 

, / HPl9l-P292j , s 

4W , 9i92(Pi9i - P292) 

V 91-92 

Ji-iJ 2 = , ,. lglg2 , ((c+l + ipi)(c-l + ipx)(2piq 1 -3q2Px-p2Q2)ql 

/ i(Pi9i -P292) , , v 

a\/ (9i - 92j 

V 91-92 

+ (2i(pi +P2)c - 4pip 2 )(pi - p 2 )qfq 2 - (c + 1 + ip 2 )(c - 1 + ip 2 )(2p 2 q 2 - 3gip 2 - Pigi)g 2 • 
The corresponding separated relations are denned by equation with the real coefficients 

a 2 

<i>{q,fi)^fi{c 2 ~l)q 3 + (2cfi 2 -b)q 2 + {^-H lf i + H 2 )q+— = 0, q = q h2 , p = i gi.api.a , (3.30) 

and quadratures in differential form look like 

9i , 92 



gi ((c 2 - l)q 2 + 4cgi/ii + 3^ 2 - i?i) q 2 ((c 2 - + 4cg 2 ^ 2 + 3/i 2 - H x 
Mi9i , M292 



= 0, 



9i 



((c 2 - l)g 2 + 4cg x /xi + 3a* 2 - ffi) g 2 ((c 2 - l)gf + 4cg 2 ^ 2 + 3/; 2 - ff x ) 



3.4.1 Deformation of the Dullin-Matveev system 

Substituting q = q\ %2 and fi — iqi, 2 pi, 2 into the non-hyperelliptic algebraic curve of genus three defined 
by the following equation 

a 2 

$1 = il(c 2 - l)g 3 + (2c M 2 + dn - b)q 2 + (/1 3 + e/1 2 - Hm + H 2 )q + — = , (3.31) 

and solving a pair of the resulting equations with respect to Hi 2 one gets deformation of the initial 
Hamilton function (|3.28[) 



Hl = Hl -Ue d —+ {ce -/ )x \ +e )j 3 + - ice - d 2 +X f (3.32) 

2 V c + x 3 x\+x 2 2 ) \(x\+x\) 



after the generalized shift (|1.8|l at 

, _ i(ce-d + a;3e) 
2^1 
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Using the same birational change (|3.13j) the curve (|3.19[) can be transformed to the canonical trigonal 
form p.6p at c 2 — c 3 = 0. 

In this case equations of motion read as 



E — = o- 

k=1 <?fc(3Mfc + 4cg fe ^ fe + 2efi k + q\(c 2 - 1) + dq k - H{) 



= iqkPk 



(3.33) 



Mfc9fc 



E 

fc=1 g/s(3/i| + 4cq k iJ k + 2e/i fe + q\{c 2 - 1) + dq k - Hi) 

3.5 Case 5 from [19] 

Let us consider last integrable system from |19j with integrals of motion 

8x3 + 5 \ 2 ax\ 



16 8(x 3 + l)J d (X3 + I) 3 / 4 V^+T 



(3.34) 



ax\J 3 



8 A 2 1 * K * ' 1 4(0:3 + I) 3 / 4 ' 
The corresponding variables of separation q± t 2 and pi j2 (|3 . 1 13.3[) are obtained from 

(3x 3 + l)J 3 i(x 2 Ji - xi J 2 ) a(a; 3 + l) 1 / 4 — iJ 3 

9l + 92 = . , r + . . r , 9l92 = ^7 . r , Pl,2 = 2lV#3 + 1 + . 

4V#3 + 1(1 — 353 J V35 3 + 1(1 - X 3 J 2(Xl + 1352j 9l,2 

Inverse transformation looks like 

a(x 3 + l) 1/4 T 19192(^1-^2) 

051 + i.£2 = „ , J3 = , 

2<?i<?2 9i - 92 

2gig2(x§ - 1) . _ ia(pi - p 2 ) a(gi + 92) 

Xl a(x 3 + I) 1 / 4 ' Jl + U2 8 (gi - 92)^3 + I) 3 / 4 2q iq2 (x 3 + 1)V* ' 

j _ ii7 = i£ig|(7^3 + l)(pi ~ Pg) 2g lg2 ( gi + g 2 )(x 3 - l)(x 3 + l) 1 / 4 
112 2a(q 1 - (?2 ((x3 + l) 1 /4 a ' 1 ' °> 

where 

_ (pi9i -P2q 2 ) 2 _ 
X3 ~ A{qi-q 2 y 

Separated relations are defined by 

a 2 i 
<S>(q,H) = ~2fiq 3 -bql + ifl 3 -H lh l + H 2 )q + — = 0, 9 = 91,2, M = « 5l,2Pl,2 (3-36) 

and we have the following quadratures in differential form 

91 92 



= 0. 



q 1 (3 f i 2 1 -H 1 -2q 2 1 ) q 2 (3i4-H 1 -2q%) 

Migi M292 I 

q 1 {^\-H 1 -2ql) 92(3/x| - Hi - 29|) 2 



(3.37) 
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3.5.1 Deformation of the system (|3.34p 

Let us add three terms to initial trigonal curve of genus three (|3.36p 



,2 



$1 = (c - 2fj)q a - (d/i + b)q z + (p* + eiif - H lf i + H 2 )q + — = . (3.38) 



The corresponding deformation of the initial Hamilton function (|3.34p has the form 



e c + 2e c d(x\ + Ax 3 + 3) 



1 



4(1 -x 3 ) V VTT^ 

after canonical transformation (jl.8l) at 

ic + ie(l + £3) id 



4 2(1 -x 3 ) 4(l + a:3) 4(l-x 3 )(l+a;3) 3 / 2 



eVl + x 3 + d+ ) , (3.39) 



The corresponding quadratures are defined by 

2 



E 



k=1 q k (3(j, 2 k + 2e(i k - 2ql - dq k - Hx) ' - 

fJ-kqk i 



E 



q k (3nl + 2e(i k - 2q 2 k - dq k ~ Hx) 2 



(3.40) 



Non canonical transformations (|3.10p relate this equations (|3.40p with similar equations (|3 . 33[) for the 
deformed Dullin-Matveev system. 



4 Conclusion 

In [2U [221 122] some new variables of separation for various integrable systems on the sphere with 
higher order integrals of motion have been obtained by brute force method. In [221 123] we introduce a 
concept of natural Poisson, which allows us to understand the geometric origin of this method and to 
find some common attributes of the variables of separation for the Kowalevski top, Chaplygin system, 
Goryachev-Chaplygin gyrostat, Goryachev and Dulllin-Matveev systems etc. 

In this more technical paper we continue our investigations in order to explicitly describe canonical 
transformations of initial physical variables to variables of separation and vice versa, to calculate the 
corresponding quadratures and to discuss possible integrable deformations of these systems associated 
with genus three hyperelliptic and non-hyperelliptic algebraic curves. 

In Section 2 we consider real variables of separation for which the separation relations have the 
real coefficients only. In Section 3 we discuss complex variables of separation and the separation 
relations with the real coefficients as above. Similar complex variables satisfying to the real separated 
equations for the Kowalevski top and Goryachev-Chaplygin gyrostat have been found in [TJ], for the 
Kowalevski- Goryachev- Chaplygin gyrostat in [TH] and for the Steklov-Lyapunov system in [25]. These 
and other known complex variables lying on the real algebraic curves are discussed in the book [2]. 

Further inquiry is related with numerical, algebro-geometric and topological analysis of the ob- 
tained quadratures. For dynamical systems associated with the (3,4) trigonal curve (|3.6p we also want 
to discuss an application of the Kowalevski-Painleve criteria to these systems, because in generic case 
solutions of the corresponding quadratures are non-meromorphic functions of time. 

We would like to thank A.V. Borisov and Yu.N. Fedorov for helpful discussions. 
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